preprint KCL-PH-TH-2012-22 



VT- Symmetric Interpretation of Double-Scaling 

Carl M Bender")!, Moshe Moshe 6 |H Sarben Sarkar c [j}] 

"Department of Physics, Washington University, St. Louis MO 63130, USA 
^Department of Physics, Technion - Israel Institute of Technology, 32000 Haifa, Israel 
c Department of Physics, King's college London, Strand, London WC2R 2LS, UK 

Abstract. The conventional double-scaling limit of an 0(A)-symmetric quartic 
quantum field theory is inconsistent because the critical coupling constant is negative. 
Thus, at the critical coupling the Lagrangian defines a quantum theory with an upside- 
down potential whose energy appears to be unbounded below. Worse yet, the integral 
representation of the partition function of the theory does not exist. It is shown that 
one can avoid these difficulties if one replaces the original theory by its 'PT-symmetric 
analog. For a zero-dimensional 0(A r )-symmetric quartic vector model the partition 
function of the PT-symmetric analog is calculated explicitly in the double-scaling 
limit. 
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1. Introduction 

The techniques of VT quantum mechanics have been used to solve several long-standing 
problems, such as the violation of unitarity in the Lee model [1] and the appearance 
of ghosts in the Pais-Uhlenbeck model [2] and in other field-theory models [3]. In this 
paper we clarify a serious problem with the double-scaling limit in quantum field theory 
that in our opinion has not yet been satisfactorily addressed; namely, that the critical 
coupling constant g cr i t is negative. Since g CT a < 0, the double-scaling limit appears to be 
unphysical because near g crit the potential is upside-down and the energy appears to be 
unbounded below. However, we show that if we approach the critical theory in a PT- 
symmetric fashion, the resulting correlated limit gives a physically acceptable quantum 
theory. Here, we study the double-scaling limit of a zero-dimensional 0(iV)-symmetric 
quartic vector field-theoretic model. (In a more detailed paper we will show that the 
same approach works for theories in higher dimension [1].) 
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To explain the mathematical problem addressed in this paper we consider the 
partition function Z(g) for a toy zero-dimensional quartic quantum field theory model: 

Z{g) = ^= tee-*'*-**'*, (1) 

V Z7T ,/_oo 



where the coupling constant g is assumed at first to be positive. The weak-coupling 
series for Z(g) is obtained by expanding the integrand of Q in powers of g and then 
integrating term by term. (This is a trivial application of Watson's Lemma [5].) The 
result is the (divergent) asymptotic-series representation for Z(g): 

n=0 

This series alternates in sign when g > and it can be Borel summed to recover the 
integral representation (Tip _ 

The integral representation ([I]) for Z(g) ceases to exist when g is negative. Thus, 
to study the behavior of Z(g) for negative g we evaluate the integral for positive g in 
terms of a parabolic cylinder function [6]: 

Z(g) = (2g)-^e 1 ^D_ 1/2 (l/^2g) . (3) 

We conclude from this formula that Z(g) is a multiple-valued function of g\ it is defined 
on a four-sheeted Riemann surface and is complex when g is negative. Thus, to evaluate 
Z(g) for g < we must define precisely the path from positive to negative g. We obtain 
four possible values for Z(g) depending on how we rotate in the complex- g plane from 
g > to g < 0. 

Even though Z(g) is a four-valued function of g, each term in the asymptotic series 
(J2j) is single valued. The resolution of this apparent discrepancy lies in identifying the 
wedge-shaped region in which a series is asymptotic to the function that it represents; 
this region is called a Stokes wedge. (For a detailed discussion of Stokes wedges see 
Refs. [5j |7j.) The angular opening of the Stokes wedge for D v (x) for x » 1 is 
|argx| < 37r/4 [5]. Thus, in the complex-^ plane the Stokes wedge is centered about the 
real-g axis and the full opening angle is 3n. This means that the series ^ is asymptotic 
to Z(g) in ^ on the negative-^ axes argg = ±7r as \g\ — > 0. Thus, we have the 
paradoxical result that on the negative axes ^ is real while Z(g) is complex. This 
happens because the imaginary part of Z(g) is subdominant (exponentially small); it is 
of order e 1//fil , as can be determined by a saddle-point calculation. 

In this paper we are interested in the partition function for negative coupling 
constant [HJ. However, we will not be interested in the complex partition function that 
one obtains by rotating x 2 + gx A from positive g to negative g in the complex-^ plane. 
Rather, we will keep the coupling constant g fixed and perform the PT-symmetric limit 
of x 2 + gx 2 (ix) 6 as e goes from to 2 [9J. The resulting partition function is 

Z(g) = ^= [ dxe~ x2 l 2+9xi /\ (4) 
c 
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where C is a comlex contour of integration. The original integration path lies on 
the real axis when e — 0. The path rotates downward into the complex plane as 
e increases. When e = 2, the final contour C comes inward from x = oo in the 
45° wedge —7n/8 < argx < — 5n/8 and goes back out to x = oo in the 45° wedge 
— 3n/8 < argx < — n/8. As a consequence, this partition function is real. This is the 
iV-dimensional version of the PT-symmetric theory that is studied in this paper. 

This paper is organized as follows. Correlated limits are explained in Sec. [2] using 
several illustrative examples. Then the uncorrelated and correlated expansions of an 
0(iV)-symmetric vector theory in zero dimensions are investigated in Sees. [3] and |4] 
and problems with implementing the double-scaling limit are described. To solve these 
problems the vector theory is reformulated as a PT-symmetric model in Sec. |5j and 
the uncorrelated and correlated expansions of this PT-symmetric model are derived in 
Sees. [6] and [7j Brief concluding remarks are given in Sec. [8j 

2. Correlated limits 

Correlated limits arise frequently in physical problems when there are two parameters, 
say e and a, and e is treated as small (e <C 1) so that it plays the role of a perturbation 
parameter. If we treat a as fixed, the solution S(e, a) to the problem is a formal 
perturbation series in powers of e: S(e,a) ~ J2^=o a n{.oi)e n . Usually, the perturbation 
series is a divergent asymptotic series in which each term in the series is negligible 
relative to the previous term as e tends to 0; that is, a n (a)e n <C a n ^.i(a)e n "~ 1 as e — > 
for all n. A correlated limit occurs when we do not treat a as fixed, but rather allow it 
to tend to a limit as e — > 0; that is, we take a to depend one: a = a(e). 

A nontrivial correlated limit arises if we choose the functional dependence so that 
all terms in the perturbation series become comparable as e — > 0. When this happens, 
the series undergoes a transmutation in which it depends on just one parameter, which 
we call 7. In this correlated limit the perturbation series still diverges, but we sum 
the series for 5(7) by using Borel summation. Correlated limits are remarkable in that 
5(7) is a universal function that reveals the essential features of the problem while being 
insensitive to specific details. Often, S(j) is entire (analytic for all 7). 

This paper examines the correlated limit of 0(iV)-symmetric quantum field theories. 
Such theories have been used to model a variety of physical phenomena because the 
large- N expansion in powers of 1/N often reveals the phase structure of the theory; 
quantities such as masses and Green's functions can be expressed as asymptotic series 
having the form J2kL a>kN~ k [10J. In the double-scaling limit iV — > 00 and g —> g CT n 
in a correlated fashion in which all terms in the 1/N expansion are of comparable size. 
In the correlated limit the sum of the series is no longer dominated by early terms; the 
N~ k power in the kth term is compensated by a sizable coefficient aj~ [H]. Furthermore, 
this limit is characterized by a universal function of the parameters that describes the 
correlated limit [11]. In the correlated limit 0(iV)-symmetric vector models represent 
discretized branched polymers. Additionally, dynamically triangulated random surfaces 
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summed on different topologies can be represented by matrix models in the double- 
scaling limit [3 [12] . 

We have devised three elementary examples to illustrate correlated limits: 
Example 1: Behavior of a nonuniformly convergent Fourier sine series near the 
boundary of its interval of convergence. Here, the number of terms iV in the partial 
sum of the Fourier series is correlated with the distance x to the boundary. This limit 
is described by the Gibbs function G(j) = Si (27) (the sine-integral function), where 
iV — > 00, x — > 0, and 7 = Nx. The Gibbs function is entire and it is universal because 
it describes the behavior of the Fourier sine series for any differentiate function f(x) 
such that f(0) 7^ and/or /(7r) 7^ 0. To explain the famous universal 18% overshoot 
exhibited by all nonuniformly convergent Fourier series at the boundary, we simply 
verify that G'(tt/2) = and that G{ir/2) = 1.18 . . .. 

Example 2: Laplace's method for the asymptotic expansion of integrals. To find 
the large- iV behavior of the Laplace integral 

POO 

Z(N)= / dre- NS{r \ (5) 



we assume that S'(r) > for all r > and use repeated integration by parts to obtain 
the complete asymptotic expansion of Z(N) as iV — > 00 [5]: 



oo 



k 


1 d ' 


fc-i 






S'{r) dr _ 


S'{r) 


r=0 



(6) 



Z(N)~ e - NS ^Y, N ~ 

k=l 

This is an uncorrelated large- iV expansion. 

Laplace's method emerges as a correlated limit of integration by parts: Suppose 
now that S'(0) is small [but that the higher derivatives of S(r) are not small at r = 0]. 
As 5"(0) 0, the kth term in ^ is approximated by 

N- k [-2S ,, (0)} k - 1 [S , (0)} 1 - 2k T(k - l/2)/r(l/2) (7) 

because this has the greatest number of powers of S"(0) in the denominator. Consider the 
correlated limit A^ — > 00, S"(0) — > 0, where ■y 2 = N[S'(0)] 2 /S"(0) is a fixed parameter. 
[We assume that S"(0) > so that 7 2 > 0.] In this limit @ becomes 

This series diverges, but we can obtain its Borel sum in terms of the parabolic cylinder 
function D„(z) [B]: 



Z(i) ~ exp ( 7 2 /4) D^)/^NS"(0). (9) 

The function Z(j) is entire. Also, it is universal because it depends only on the two 
numbers 5(0) and S"(0), and thus it applies universally to all functions S(r) with these 
particular values. [In contrast, the uncorrelated series ^ depends on all derivatives of 
S(r) at r = 0. ] 

For the special value 7 = 0, D_i(0) = a/vt/2 gives the famous result 

Z(N) ~ e- JVS(0 V 7r /[2WS"(0)] (N — >• 00) (10) 
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of Laplace's method applied to (pH). The asymptotic formula (10) is a limiting case of 
the correlated limit ^ for which 5"(0) = and S"(0) > 0. Thus, ^ describes the 
approach of Z(j) to Laplace's asymptotic formula (10). To be precise, the asymptotic 
expansion of the Laplace integral is controlled by the function S(r) at its Laplace point 
(that is, its maximum point). Laplace's formula (10) gives the asymptotic behavior 
of the integral for the special case in which S(r) is level (has a vanishing derivative) 
at its Laplace point. Thus, the correlated limit describes in a smooth and universal 
fashion what happens as the derivative of S(r) approaches at the Laplace point, just 
as the Gibbs function describes in a smooth and universal fashion how a nonuniformly 
convergent Fourier series for f(x) behaves as x approaches the boundary of the interval. 

Example 3: Transition in a quantum-mechanical wave function between a 
classically allowed region and a classically forbidden region. This transition is described 
by the universal Airy function Ai(7), which is obtained by summing the divergent WKB 
series in the correlated limit h — » 0, x — > (where x is the distance to the turning point), 
with the ratio 7 = x 3 ^ 2 /h held fixed. Consider the one-turning-point problem for the 
Schrodinger equation h 2 <p"(x) = Q(x)<f)(x), where Q(x) is linear in x near the turning 
point at x 



0: Q(x) ~ ax (x — >■ 0). When x 7^ the WKB approximation to (f)(x) is 

(fr^O), 



Vkb [ x ) 



cxp 



where S (x) = ±^/Q{ 
relation 2S (x)S,, 



x 



SAx) 
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obeys the recursion 



x 



-SL 



n-l 



X 



n=0 

= -Q'(x)/[4Q(x)], and S n (x) 
Y%Zi S j{ x ) S n-j{x) for n > 2. 
The WKB approximation 4>wkb(x) is invalid at the turning point because it blows 
up at x = while the exact solution <p(x) remains finite. To see how the WKB solution 
behaves as x approaches 0, we choose the negative solution for So: So(x) ~ —y/ax. 
Then, Si(x) ~ — l/(4x) and the functions S n (x) have the asymptotic form 



S n 



n 1/2-71/2^1/2-371/2 



X 



2 s n (x 0), 



(12) 

where s n = (4 — 3n)s n _i — \ YTj=\ s j s n-j for n > 2 and s = Si = 1. The numerical 
coefficients s n become complicated as n increases, but an easy way to understand the 
WKB series is to exponentiate it: 



exp 



00 

n=2 



.n—1 



dx S xi 



x 



fc=0 



-i)W 



(13) 



where 7 = a 1 / 2 i"^/ 2 
now have the simple form g^ 



(!''"x"'"/h. The numerical coefficients g n in the exponentiated series (13) 



7i- 1/2 r(3A; + l/2)9" fc /r(2fc + 1). The series on the right 



side of (13) is just the asymptotic expansion of the Airy function Ai(z) as z — > +00: 



Ai(z) 



-1/4 



2A 



e 3 Z ^{-lfg k z 



-3k/2 



(14) 



fc=0 



This series diverges but it is Borel summable; in the correlated limit h — > 0, x — > 0, 7 
fixed, the solution to the Schrodinger equation is 0(7) = cAi(7), where c is an arbitrary 
constant. This shows that the solution to the famous one-turning-point problem is a 
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correlated limit. The Airy function is an entire function of 7 and it is universal because 
it is valid for all potentials Q(x) that vanish linearly at the turning point. 



3. Uncorrelated large-iV series for quantum field theory in zero dimensions 



The partition function 



d N+1 



xexp 



N+l 



n=l 




(15) 



represents a zero- dimensional quartic quantum field theory having 0(N + 1) symmetry. 
The coupling constant A is assumed to be positive so that the integral converges. To 
derive an uncorrelated large-iV expansion of Z, we exploit the rotational symmetry by 
introducing the radial variable r, Y2n=i x n = Nr 2 , and we let A = g/N. The partition 
function now takes the one-dimensional form 



A 



N+l 



dre~ NL ^\ 



(16) 



where A N = 2tc n ^ 2 /T(N/2) is the surface area of an iV-dimensional sphere of radius 1 
and L(r) = r 2 /2 + gr 4 /4 — logr. We emphasize that we must assume that g is positive 



so that the integral representation (16) for the partition function Z converges. 



The integral ( 16 ) is a Laplace integral and the correlated expansion of the Laplace 
integral ^ was considered earlier. However, in contrast with Z(N) in ([5]), L'(r) is not 
positive. We will see that as a result, the correlated limit of (flBl lies in a different 



universality class and is characterized by a different universal function; specifically, the 



universal function for (16) is an Airy function rather than a parabolic cylinder function. 



The standard procedure (Laplace's method) for finding the large-iV asymptotic 



behavior of the integral (16) begins by locating the Laplace points, which are the zeros of 
L'(s) = r + gr 3 — 1/r. Just one Laplace point r = y(G — l)/(2g), where G = y/1 + Ag, 
lies in the range of integration < r < 00. Laplace's method relies on the crucial 
fact that if ro is a global minimum of L(r), then as N — > 00 the entire asymptotic 
expansion of the integral is determined by the behavior of L(r) on the infinitesimal 
region r — e < r < r + e containing r . [The Laplace point r is a global minimum 
because L"(r) = 1 + 3gr 2 + 1/r 2 and L" (r ) = 2G > 0.] We thus Taylor expand L(r) 
about r and conclude that for large N, 



ro+e 



Z~A N+l e~ NL ^ I dre- NG ^ 



exp 



ro-e 



N 



E M L(fc) (ro) (r - ro) ' 



fc=3 



To evaluate this integral we translate the integration region by substituting t = 
r — ro and then we perform the scaling t = u/y/ NG. Because the sum in square 
brackets begins at k — 3, all terms in the sum are small as iV — > 00. Thus, we can 
expand the exponential in powers of u. Making only transcendentally small errors, we 
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extend the range of integration to — oo < u < oo, and perform the Gaussian integrals 
term-by-term. This yields an infinite series in powers of 1/N: 



A N+1 e 



-NL(r ) 



7T 



(17) 



fc=0 



This is the full asymptotic expansion of Z and it is the uncorrelated large- N expansion 



of the partition function (15). Only integer powers of iV 1 appear in this expansion; 



there are no half-odd-integer powers because Gaussian integrals over odd powers of u 
vanish. The first three coefficients in this series are explicitly 

5 - 6G 2 - G 3 385 - 924G 2 - 10G 3 + 684G 4 + 12G 5 - 143G 6 



a 



2AG 3 



a>2 



1152G 6 



4. Correlated limit of the large- A" expansion (17) 



Proceeding formally, we now attempt to construct a correlated large-iV limit of the 



expansion in (17). For all terms in the expansion to have the same order of magnitude, 
the correlated limit must be N — > oo and g — > g crit = —1/4 (that is, G — > 0) with 
7 = NG 3 /2. In this limit the asymptotic approximation (17) becomes 



-NL(r ) 



^NG/tt 



1 + 



+ 



385 



487 46087* 



+ . . . 



We recognize that the series in (18) is the asymptotic expansion for large 7 of 



Bi(7 2 / 3 ) v /7re 2 ^/ 3 7 1 / 6 . [This is like the series in (13) for the Airy function Ai, but it 
does not alternate in sign.] Thus, we are tempted to conclude that the correlated limit 



of the series (17) is 



Z ~ A N+1 e NL ^2- 1 /%N- 1 ^Bi( 1 2 ^)e- 2 ''/ 3 . (19) 
However, this correlated limit is invalid because it requires that g < 0. Furthermore, 



to obtain the Airy function Bi in (19), we have had to sum a nonalternating divergent 



series; such a series is not Borel summable. 



5. PT-symmetric reformulation of the theory 



In place of (15) we consider the 0(A^ + l)-symmetric partition function Z = 
Re J d N+1 x e~ L , where we assume that is an even integer. The Lagrangian L has 
the form 

3=1 v i=i 7 

The multiple integral above is taken on the real axis and it converges if e < 1. 

We let A = gN~ £ l 2 and again introduce the radial variable r by ^2^=i x n = Nr 2 . 
The crucial assumption that A^ is even allows us to extend the radial integral to the 
entire real-r axis: 

1 r°° 

Z = -A N+1 dre- N£ < r \ (21) 
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Figure 1. The four saddle points [zeros of C'(r)] for the case g < 1/4 with the up 
and down directions indicated by the letters U and D. The steepest-descent path is 
■pT-symmetric (left-right symmetric) and runs along the real axis between the distant 
pair of saddle points. Then it turns downward, curves, and asymptotes at — 7r/4 and 
at — 37r/4. For large N the contribution to the integral is localized at the pair of saddle 
points nearest the origin. 



where C = r 2 /2 + gr 2 (ir) £ /(2 + e) — logr. This Laplace integral is real because the 
integrand is VT symmetric; that is, the integrand is invariant under r — > —r and 
i — >• —i |9J. Without its logarithm term, £ has a standard "PT-symmetric structure that 
has been studied in great detail [9]. The logarithm in the exponent does not make the 
integral complex. [We can make C explicitly VT symmetric by including an additive 
constant logr — > log(zr) and taking the branch cut to lie on the negative- imaginary axis. 
The additive constant log* has no effect on the forthcoming steepest-descent analysis.] 



To achieve the PT-symmetric partition function (21) we have had to work in a space of 
odd dimension N + l. This requirement can be traced to the absence of a distinct parity 
operator in even-dimensional space, where the sign of x can be changed by a rotation. 
To obtain a quartic theory as e — > 2, we must redefine the boundary conditions on 



the integral (21) accordingly: For any e > the integral converges if the integration 
contour lies inside a pair of PT-symmetric Stokes wedges centered about — tie/ (4 + 2e) 
and — (47r + 7re)/(4 + 2e). The wedges have angular opening vr/(2 + e) and contain 
the real-r axis if e < 1. As e increases above 1, the wedges rotate downward into the 
complex plane and become narrower. At e = 2 the wedges are centered about — 7r/4 
and — 37r/4 and have angular opening 7r/4. 



6. Uncorrelated limit of the PT-symmetric theory (21) 



To find the large- A" asymptotic behavior of Z in (21) when e = 2, we use the method 
of steepest descents j5]. The saddle points are the zeros of C(r) = r — gr 3 — 1/r. There 
are four saddle points, which are the roots of r% = (1 ± y/1 — Ag)/(2g). If g < 1/4, the 
saddle points are all real and are shown in Fig. [TJ (If g > 1/4, the saddle points are all 
imaginary, but this does not affect the asymptotic analysis of the double-scaling limit.) 

We determine the directions of the saddle points by calculating C"{r) = 1 — 3gr 2 + 
1/r 2 . At the saddle points C" = ^=2^/1 — Ag. For the minus (plus) sign the steepest- 
descent path moves away from the saddle point in the imaginary (real) direction. Thus, 
the complete steepest-descent path, as shown in Fig. [TJ follows the real-r axis until it 
reaches the distant pair of saddle points. Then it turns downward and curves off at the 
angles —ir/A and — 3n/4. The entire contribution to the uncorrelated limit comes from 
an infinitesimal region surrounding the two saddle points closest to r = 0. 
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7. Correlated limit of the "PT-symmetric theory 



In the correlated limit the critical point is determined by requiring that C"(r) = in 
addition to C'(r) = 0. The critical value of the coupling constant is g CT ft = 1/4. At this 
value of g there is a coalescence of the two uncorrelated saddle points at r crit = ±\/2, 
as shown in Fig. |2j The steepest curve is shown in Fig. [2] and the contribution to the 
integral in the double-scaling limit comes from four infinitesimal regions at the saddle 
points, two horizontal and two at 60° angles. We must now evaluate four convergent 



integrals of the form J Q dr e° 



] -br i 



where L = oo or L = ooe 17r ^ 3 . Each of these 



integrals separately gives a combination of Airy functions and a 2-F2 hypergeometric 
function. However, we will now demonstrate that when the four integrals are combined 



there is a dramatic simplification and the result is identical to that in (18). 




Figure 2. Correlated limit of Fig. [T] At g = g cr i t the quadratic saddle points coalesce 
to form cubic saddle points with up and down directions indicated by U and D. The 
steepest path now runs along the real axis between the saddle points. At the saddle 
points it veers off at 60° angles, curves, and asymptotes at — 7r/4 and at — 37r/4. 

We begin with the integral f dre~ r2+ri '/(Vt)^ where 7 = G 3 N/2 and G = y/1 — Ag. 



Note that the sign of A is reversed in (21) compared with (15). There are four saddle 
points of cubic type. This gives rise to two integrals for each saddle point: 

dr + dr + complex conjugate. 

Jo 

In the first integral we replace r — > — r and in the second we replace r — > re~ l7T ^ 3 . The 
first integral then becomes 3 1 / 3 7 1 / 6 J °° e^" 2- " 3 , where / = — 3 2 / 3 7 1//3 , and the second 
becomes 3 1 / 3 ^ 1 / 6 e - :(7r / 3 e fu 2 -u 3 ^ w h ere j = e i7r/3 32/3^1/3 _ N ex ^ we use the identity 

2 4 4/ 3 ^ 



3' 3' 27 



where the generalized hypergeometric function is defined by the Taylor series 



(22) 



r («) r ( & ) ^ ™ !r ( c + n ) T ( d + n) 



For both the first integral and the second integral, z = —A'j/3. Thus, the two 
hypergeometric functions are real and identical, their sum exactly cancels, and only the 
Bi functions remain. There are four such Airy functions: 

27re- 2 ^ 3 3- 4 / 3 [2Bi ( 7 2/3 ) - e 2 -/ 3 Bi ( e 2 -/ 3 7 2 / 3 ) - e - 2l7r / 3 Bi ( e - 2 -/ 3 7 2 / 3 )] . 



VT -Symmetric Interpretation of Double- Scaling 



10 



We simplify this expression by using the identity Bi(z) + cuBi (uz) + u 2 Bi (uj 2 z) = [6]. 
The final result for the integration is therefore 27re -27/ ' 3 3~ 1 / 3 Bi (y 2 ^ 3 ). This concludes 
the demonstration. 

8. Final remarks 

We have shown that the key to constructing and interpreting the double-scaling limit 
is continuing in e rather than in g. In the past the way to deal with a negative-quartic 
potential has been to imagine an analytic continuation from positive g to negative g. 
However, this procedure would give complex-energy eigenvalues because the path from 
positive to negative g crosses branch cuts in the coupling-constant plane, which is a 
multisheeted Riemann surface [13]. We have used a different analytic continuation, 
e : 0—7-2, which does not involve the coupling constant. To do so, we have written 
the potential in manifestly PT-symmetric form. The complex potentials for all e > 
give real partition functions and at e = 2 we have made contact with the quartic theory 
obtained in the conventional double-scaling limit. 
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